Convective instability in inhomogeneous media: impulse response in the subcritical 

cylinder wake 
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We study experimentally the impulse response of a cylinder wake below the critical Reynolds 
number of the Benard-von Karman instability. In this subcritical regime, a localized inhomogeneous 
region of convective instability exists which causes initial perturbations to be transiently amplified. 
The aim of this work is to quantify the evolution resulting from this convective instability using two- 
dimensional particle image velocimetry in a hydrodynamic tunnel experiment. The velocity fields 
allow us to describe the evolution of wave packets in terms of two control parameters: the Reynolds 
number and the magnitude of the imposed perturbation. The temporal evolution of energy exhibits 
a transient algebraic growth at short times followed by an exponential decay. 
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I. INTRODUCTION 

Flow past a circular cylinder is a classic prototype for 
studying hydrodynamic instabilities and bifurcations in 
separated flows pHZI]. Moreover, the academic case of 
a two-dimensional cylinder wake can be used as a basic 
model for many real situations, including the flow behind 
support cables or around an airfoil. When the Reynolds 
number Re = UqD/p (where Uo is the free-stream veloc- 
ity, D is the cylinder diameter, and v is the kinematic 
viscosity) reaches a particular critical value (Re c » 47 
for an infinitely long cylinder, [H, 0), a sustained pe- 
riodic shedding of opposite-signed vortices gives rise to 
the well-known Benard-von Karman vortex street. 

The transition to sustained oscillations can be de- 
scribed locally, via wake profiles at different spatial sta- 
tions, or globally, viewing the 2D wake as whole. To 
elaborate further, we first recall the standard distinction 
between convective and absolute instability in parallel 
flows illustrated in Figs.QJa) andQJb). An instability is 
convective if a perturbation grows but is simultaneously 
advected with the flow such that the disturbance decays 
at any fixed point, as in Fig. [Ha), while it is absolute if 
the perturbation grows at a fixed spatial position, as in 
Fig- H^b) [cl-Hoj. A useful way to distinguish convective 
from absolute instability is in terms of propagating edge 
or front velocities: assuming a positive leading-edge ve- 
locity V + , as is the case here, convective instability corre- 
sponds to a positive trailing-edge velocity V~ , while ab- 
solute instability corresponds to a negative trailing-edge 
velocity 0-GJ]. 

Local analysis of wake profiles, i.e. a parallel flow ap- 
proximation, gives a picture of the transition to sustained 
oscillations as follows [1, [rJ [l4[ • Above Re ~ 5 there is 
a local region of convectively unstable flow in the wake. 
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FIG. 1: Sketch of the space-time response of flows to infinites- 
imal perturbations, (a) and (b) correspond to local analysis, 
i.e. parallel flows, illustrating the distinction between (a) con- 
vective and (b) absolute instabilities. The leading-edge veloc- 
ity V + is positive in both cases. For convective instability, 
the trailing-edge velocity V~ is also positive, while for abso- 
lute instability it is negative, (c) and (d) correspond to the 
global wake illustrating the distinction between the (c) sub- 
critical, Re < Re c , and (d) supercritical, Re > Re c , cases. In 
the subcritical case, the perturbation reaches a maximum and 
subsequently decays. In the supercritical case the perturba- 
tion continues to grow until it saturates nonlinearly. The edge 
velocities can be obtained over finite times only in experiment 
(see Sec. UTI}. 



Above Re ~ 25 there is additionally a pocket of abso- 
lutely unstable flow within the region of convective in- 
stability. Once the pocket of absolutely unstable flow 
becomes sufficiently large, the flow becomes globally un- 
stable. Globally, however, one does not observe the onset 
of a locally absolutely unstable region. Instead, one finds 
that below Re c the wake responds to perturbations, but 
only transiently, as perturbations are advected through 
the systems. This is illustrated in Fig. [He). Above 
Re c , as in Fig. [T|d), perturbations grow and lead to a 
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FIG. 2: (a) Schematic diagram of the experimental setup in 
the hydrodynamic tunnel, (b) Sketch showing a top view of 
the cylinder during an impulsive rotation of intensity U p . 



synchronized wake in the formation of what is called a 
global mode 0, HH . Hence, even though convective and 
absolute instability are strictly defined for parallel flows 
and streamwise period flows 0, [ToL [f6j , inhomogeneous 
flows, such as flow past a cylinder, may exhibit similar 
characteristics. The transient response in the subcritical 
regime has the hallmarks of convective instability while 
the global instability above Re c has the hallmarks of ab- 
solute instability. 

Despite the large body of work on the cylinder wake, 
only a few experimental studies have examined the wake's 
subcritical behavior [l7| and none have quantitatively 
characterized the transient dynamics in this regime. The 
goal of the present paper is therefore to analyze quan- 
titatively the subcritical regime using a well-controlled 
experiment. 

After describing the experimental setup, we focus 
firstly on the evolution of amplified wave packets and 
on obtaining the leading-edge, trailing-edge, and group 
velocities in an experimental setting. The decrease of 
trailing-edge velocity towards zero when approaching the 
global instability threshold confirms the transition to an 
absolute instability. The subcritical behavior is further 
characterized in terms of the evolution of the maximum 
amplitude of the wave packet, and its space-time posi- 
tion, as a function of the strength of the perturbation 
and the distance to the Benard-von Karman instability 
threshold. Finally, to quantify the transient growth phe- 
nomenon due to the inhomogeneity of the media, we an- 
alyze the temporal evolution of the energy of the pertur- 
bation. 



II. EXPERIMENTAL SETUP 

A cylinder of diameter D = 5 mm is placed in a hydro- 
dynamic tunnel of section 100mm x 100mm (see Fig. [5]) 
with a nearly plug flow in the test section. The bound- 
ary layer width in the tunnel walls is of ~ 10 mm in the 
region of interest. The cylinder span is 98 mm which 
covers practically the whole height of the tunnel. We 
use a Cartesian coordinate system, placed in the cylin- 
der center, with the x-axis pointing in the flow direction 
and the z-axis running along the cylinder centerline. The 
cylinder can be put into rotation to provide impulse per- 
turbations to the flow. Measurements are taken in the 
horizontal mid-plane of the channel. 

We define and work with the reduced Reynolds number 
e = (Re — Re c )/Re c , corresponding to the distance from 
the global instability threshold. The critical Reynolds 
number measured in this experiment is Re c w 64, which 
is larger than the ideal 2D case, primarily due to con- 
finement and blockage effects. We adjust the Reynolds 
number in the experiment by controlling the flow rate in 
the tunnel. In the results that follow, e varies from —0.30 
to —0.04. For this range of Reynolds numbers, the flow 
remains two-dimensional over nearly the entire cylinder 
span. 

Impulse perturbations consist of applying very short 
rotary motions to the cylinder. Rotation is controlled by 
a programmable microstepping electronic module which 
gives a resolution of 1/256 per full step, allowing for a 
precise control of the cylinder motion. In the present 
work, we fix the non-dimensional time interval over which 
the cylinder is rotated, and use the speed of rotation 
as the amplitude of the perturbation. More specifically, 
the perturbation amplitude is defined from the tangen- 
tial speed of rotation at the cylinder surface U p (see 
Fig. [2]). We consider three perturbation amplitudes given 
by three values of the non-dimensional rotation speed: 
Up/Uo = 75, 100 and 125, which we refer to as small, 
medium and strong, respectively. Uq is the measured 
velocity of the flat profile. The small-amplitude pertur- 
bation is the smallest perturbation that produced an ob- 
servable response in the wake. For all applied perturba- 
tions, the cylinder is rotated for a fix dimensionless du- 
ration given by At/T ac i v ~ 0.2, where At is the dimen- 
sional duration of the perturbation and T a d v — D/Uq. 
Note that since Uq varies with Reynolds number, both 
the dimensional tangential speed U p and the dimensional 
duration At, vary with Reynolds number. 

The wave packet generated by the perturbation can 
be observed qualitatively by visualizing streaklines, as in 
the work of [17 ;]. Fig. [3] shows this type of visualization 
compared with snapshots of the cross-stream component 
of the velocity field obtained by 2D Particle Image Ve- 
locimetry (PIV) in the horizontal mid-plane. This high- 
lights a significant difference between the present exper- 
iment and previous studies of the subcritical wake re- 
sponse. Velocity field measurements obtained by PIV 
permit one to study directly the spatial structure of ve- 
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FIG. 3: Visualization of the impulse response at three succes- 
sive time instants (from bottom to top) in the cylinder wake. 
Left: streaklines obtained from Fluorescein dye visualization. 
Right: instantaneous cross stream velocity field obtained from 
PIV measurements. 



locity perturbations as they evolve. The velocity per- 
turbations exhibit a well-defined maximum in the wake. 
These are impossible to determine from streakline records 
since the deformation amplitude of injected dye never 
decays downstream, an artifact caused by mass conser- 
vation of the dyed fluid. Therefore, streakline deforma- 
tion does not give information about the amplitude of 
the velocity fluctuations U y , and thus about the spa- 
tial evolution of the impulse response. PIV acquisition 
and post-processing have been performed using a LaVi- 
sion system with an ImagerPro 1600 x 1200 CCD camera 
with a 12- bit dynamic range capable of recording double- 
frame pairs of images up to 15 Hz and a two rod Nd:YAG 
(15mJ) pulsed laser. The time lapse between two frames 
is set to 20 ms. Finally, additional post-processing and 
analysis has been carried out with Matlab and the PIV- 
Mat Toolbox. 



III. EVOLUTION OF THE WAVE PACKET 

The convective nature of the impulsively perturbed 
cylinder wake is clearly illustrated in the experimental 
data shown in Figs. @|a) and0Jb). In Fig. Ufa), cross- 
stream velocity U y profiles, measured at successive times 
along the wake symmetry axis y = 0, are stacked up to 
form a spatio-temporal diagram. The cross-stream veloc- 
ity on the symmetry axis is the ideal quantity to use to in- 
vestigate the perturbed flow field because it is everywhere 
zero for the unperturbed flow. Figure EJb) is similar to 
Fig. 0Ja) except that the envelope [obtained through a 
Hilbcrt transform of U y as illustrated in Fig.|4jc)] is plot- 
ted. 

We measure a number of velocities associated with 
the space-time evolution of the wave packets. From the 
Hilbert transform of the cross-stream velocity we are able 



to extract, at each time instant, the envelope of the wave 
packet as illustrated in Fig. SJc) . Procedurally, we define 
the leading and trailing positions as the inflection points 
of hyperbolic tangent fits to the relevant parts of the 
envelope at each time. Such fits, together with the asso- 
ciated inflection points, are included in Fig. HJc). Then 
we obtain the leading-edge velocity V + and the trailing- 
edge velocity V~ from the collection of fronts extracted 
from the envelopes at short times (< 50T a( fo). Likewise, 
we obtain the envelope maximum at each time and from 
these data we define the group velocity V g to be the ve- 
locity of envelope maximum, as this gives the speed of the 
packet as a whole. These velocities are the inverse of the 
slopes shown in the spatiotemporal diagram of Fig. 0Jb). 
Finally, we define the phase velocity V p as the translation 
speed of the vortices in the evolving wave packet as seen 
in Fig. QJa). Note that all vortices move at essentially 
identical speed. Moreover, the leading-edge velocity V + , 
which is effectively determined by the velocity of the first 
advected vortex released from the cylinder, is the same 
as the phase velocity to within experimental precision: 
V+~V P . 

We show now that the edge velocities obtained proce- 
durally through inflections points of the envelop corre- 
spond very closely to the edge velocities defined in con- 
vectively unstable systems. (See in particular Ref. IT3 .) 
In this case one defines V~ and V + to be the x/t rays 
that separate regions of growth from regions of decay 
in the space-time diagram. The growth rate along rays 
defining the edges is zero. For the experimental data, the 
growth rate along rays is shown in Fig. HJd) , where the 
envelope of U y , sampled along different x/t rays, is plot- 
ted against ray velocity. The edge velocities V + and V~ 
obtained from inflection points are indicated by vertical 
dashed lines in Fig. [31(d) . It can be seen that these are 
in excellent agreement with the zero growth-rate rays. 
Moreover, the group velocity V g (vertical solid line) is 
also agrees very well with the maximum of the envelop 
in the ray representation. 

Thus it can be seen that the velocities are all quite well 
defined experimentally, as least up to the time at which 
the perturbation reaches its maximum, even though the 
flow is spatially inhomogcncous and the wave packet lives 
only a finite time. As expected none of the velocities 
(including V p ) are strictly constant over this space-time 
region, but they are very nearly so outside of the near- 
wake region. We have not attempted to extract their 
detailed variation in the present study. 

We now consider the behavior of the different velocities 
as |e| — > (i.e when Re — > Re c ). In Fig. |5fa) we can see 
that the phase velocity V p , and equivalently the leading 
edge velocity V + , normalized by the inflow speed, ex- 
hibit a slight decrease with increasing Reynolds number. 
Generally, the value of V p /Uq in the subcritical regime 
(e < 0) compares well with the value of 0.88 reported in 
the literature [r| for the supercritical (e > 0) regime. 
The group and trailing edge velocities are summarized 
in Fig. H^b). Both velocities decrease with increasing 
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FIG. 4: Evolution of the wave packet seen in spatiotemporal diagrams for e — —0.137. (a) Cross-stream velocity U y (x,i), (b) 
envelope of U y (x,t), (c) U y (x) and its envelope (bold dashed curve) at a fixed time (t = 24T a d v )- Also shown (thin curves) are 
tanh fits to the relevant part of the envelope. The infection points and envelope maximum are indicated, (d) Envelope of U y 
sampled along various rays as a function of ray velocity. The edge and group velocities, obtained from the inflection points and 
envelope maxima, are indicated with vertical lines. 



Reynolds number as one would expect in approaching 
the global wake instability. To within experimental ac- 
curacy, V~ approaches zero as e — > 0. 

From the spatiotemporal diagrams we are also able to 
pinpoint the position in space and time (X max ,T max ) at 
which the subcritical response reaches maximum ampli- 
tude. See Fig. EJa). Again we use the cross-stream ve- 
locity component U y . We let A max denote the maximum 
of the response, so that A max = U y (X max ,T max ). Fig- 
ures [S^b) through [IJd) show the dependence of the max- 
imum on reduced Reynolds number. It is not surprising 
that A max grows when the global instability threshold is 
approached (|e| — > 0), as can be seen in Fig.[Bfb), because 
the susceptibility of the flow increases near the onset. On 
the contrary, the behavior of X max (which diminishes 
when |e| — > 0, see Fig.^c), deserves a further comment, 
since it seems to be opposite to the case of other insta- 
bilities in the subcritical regime where the characteristic 



length scale increases when |e| — ► (e.g. the penetra- 
tion length in pre-transitional Rayleigh-Benard convec- 
tion PJI). The fact that the maximum of the instability 
moves closer to the cylinder (i.e. that X max diminishes) 
when approaching the threshold while the time at which 
this occurs remains constant [on a nondimensional scale 
normalized by the advective time scale T a d v , Fig. Efd)] is 
actually consistent with V g decreasing as the convectively 
unstable system tends to the absolute instability thresh- 
old. The difference with the case of Rayleigh-Benard con- 
vection comes from the effect of the mean flow advection, 
which modifies the physical meaning of the penetration 
lengthscale. 

For a given value of e, an increase in the strength of the 
perturbation produces a response in which the maximum 
amplitude A max is slightly larger and occurs at a position 
closer to the cylinder (smaller X max ). Approaching the 
threshold reduces the effect of changing the perturbation 
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is only transient. Such behavior is known in inhomo- 
geneous flows [H. [2p| - [27| . In should be noted, however, 
that inflow or other noise may modify the picture in that 
the localized region of instability acts as an amplifier and 
sustained dynamics may arise in some cases even in the 
subcritical regime [U M, HE IM H| • 

An increasingly common approach to quantifying the 
transient response of flows is in terms of their transient 
energy growth 28]. Such an analysis provides a global 
measure of the response dynamics. Here we undertake 
such an analysis of experimental data. 

We define the perturbation energy from our measure- 
ments as follows: 



E(t) 



dxdy 




where 



-0.3 -0.2 -0.1 

e = (Re - Re c )/Re c 



FIG. 5: (a) Phase velocity V p (which is also V + ) and (b) 
group velocity V g and trailing-edge velocity V~ as a function 
of e. The dashed line in (a) corresponds to the value of V p /Uo 
for the Benard-von Karman wake (Re > Re c ). The dashed 
lines in (b) are linear fits to the V g and V~ data. Experimen- 
tal points are labeled as •: small perturbation, o: medium 
perturbation, □: strong perturbation. Horizontal error bars 
indicate the uncertainty in the Reynolds number measure- 
ment. 



strength. 



IV. TRANSIENT GROWTH 



As noted in the introduction, an important feature of 
the cylinder wake is the inhomogeneous nature of the 
flow. Hence, even though the subcritical response just 
presented has much of the general character of convective 
instability, perturbations do not grow indefinitely (even 
linearly) as they would for a homogeneous convectively 
unstable system. Rather, a localized initial perturbation 
grows at first, due to local flow features near the cylin- 
der, but is simultaneously advected downstream into a 
region of stability where the perturbation decays. Hence 
in the absence of any inflow noise the impulsive response 



U x - U xb 



Uy ~ U, 



ybase 



where the base flow (U x base,U y base) is the measured 
steady flow before any perturbation. 

Experimentally, the total energy of the perturbation 
can only be measured while the packet is contained in the 
observation window. Hence, the energy calculated from 
the velocimetry data does not include the contribution 
from vortices that have been advected out of the mea- 
surement area. In order to quantify this effect we com- 
pare the energy calculated using two different streamwise 
sizes for the integration area A: the total energy E to t, 
where the whole measurement window is used, and the 
energy of the first half of the wave packet Ehaif, where 
the downstream boundary of the integration area A is 
given by the time-dependent position of the maximum 
perturbation amplitude X max (t) - see Fig.[7[a). We note 
that X max (t) is time-dependent and should not be con- 
fused with X max shown in Fig (He) ■ The quantities are 
related via: X max = max t X max (t) . We find that the ra- 
tio of Etot/Ehaif remains approximately constant (« 2) 
in time, which means that the measurement window is 
sufficiently large to capture the dynamics of the pertur- 
bation growth and decay before the effect of the flow 
structures advected away from the downstream bound- 
ary of the measurement window becomes significant. 

The time evolution of the energy is shown in Fig.[7jb). 
The value of E tot , as well as E ha i f and E tot - E ha i f , are 
shown. All energies are all normalized by Eq, the value of 
E t ot at the first measured instance following the impulse. 
The perturbation energy initially undergoes growth un- 
til t ss 50T a dv, at which point the energy decays. This 
is precisely the transient growth dynamics expected of 
convective instabilities in inhomogeneous media. 

The late time behavior of the energy corresponds to 
the exponential decay of the least stable normal mode in 
a stable region, e.g. @,[2!|. As the system approaches the 
absolute instability threshold (e — » 0) this decay becomes 
slower. This can be be seen in Fig. [7£c), where the decay 
rate is given by the slope of the curves of log E to t and 
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FIG. 6: (a) Definition of the maximum of the perturbation on the spatiotemporal diagram, (b) Evolution of the maximum of 
the perturbation A ma x, (c) X ma x and (d) Tmax as a function of the reduced Reynolds Number e. A max , Xmax and Tmax are 
rendered nondimensionalized respectively by the free stream velocity Uo, the cylinder diameter D and the advective timescale 
T a dv = D/Uq- Experimental points are labeled as •: low perturbation, o: medium perturbation, □: strong perturbation. 



more quantitatively in Fig. [SJa) where the asymptotic 
decay rate is plotted as a function of reduced Reynolds 
number. For comparison, decay rates from linear sta- 
bility computations Q shown with a solid curve. The 
agreement is excellent. 

The energy growth at short times can be examined 
using the maximum of the energy E max and the time 
for which this maximum is reached t max as a function 
of e. See Fig. [SJb). Consistently with the approach of 
the absolute instability threshold, E max increases when 
|e| — > 0. The time t max where the maximum energy is 
reached remains approximately constant. 



V. DISCUSSION AND CONCLUSIONS 

The convectively unstable wave packets produced by 
an impulsive perturbation in a subcritical cylinder wake 
have been studied experimentally. Velocity field mea- 



surements obtained by PIV have permitted us to charac- 
terize quantitatively the instability wave, shedding light 
on points that remained not clearly analyzed in the lit- 
erature. Firstly, probing the perturbation of the ve- 
locity field due to the instability shows unambiguously 
that there is a well-defined maximum of the perturba- 
tion amplitude in the wake. This differs significantly from 
the picture given by the streakline visualizations usually 
used to illustrate the convective instability (l7| . which 
distort the observation of the actual growth and decay. 
The value of this maximum and its position downstream 
the cylinder depend on the Reynolds number and, less 
markedly, on the strength of the perturbation, whereas 
its position in time remains constant. In addition, the 
evolution of the wave packet has been characterized with 
respect to the Reynolds number using the leading and 
trailing fronts as well as a typical group velocity. The 
measured velocities are consistent with the transition 
from a convective to an absolute global instability as 
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FIG. 7: Transient energy growth in experiment, (a) Illustration of the integration area for the energy on the instantaneous 
velocity field U y . Etot and Ehaif are calculated over the white and black rectangles, respectively. The downstream boundary of 
the black rectangle moves at a speed V g and tracks the maximum of the perturbation envelope as indicated in the upper plot, 
(b) Energy as a function of time for e = —0.14. The dots are the experimental data of Etot- The solid, dashed and dotted lines 
are smooth fits of E to t, Ehaif and Etot — E^aif, respectively, (c) Time evolution of the energy on logarithmic scale, for three 
different values of the reduced Reynolds number e: • e = —0.19, + e — —0.14 and □ e — —0.1. The slope of the energy decay 
on the logarithmic scale gives the decay rate a. 



the Reynolds number increases towards the Benard-von 
Karman instability threshold (see Fig. [3]). 

The analysis of the transient energy growth associated 
to the instability also deserves a final comment. The 
qualitative features of the temporal evolution of the en- 
ergy agree with the standard picture of transient growth 
due to convective instability in inhomogeneous media, i.e. 
a short-time algebraic growth followed by an exponen- 
tial decay at late times. However, the measured values 
of max(E to t / Eq) are remarkably low (always less than 
10) when compared to the values obtained by numerical 
computations (of order 10 3 , see e.g. [H, H3)- The main 
reasons for this discrepancy are most likely the fact that 
in the experiment one does impose an optimal pertur- 



bation and also the fact that extracting Eq, the initial 
energy of the perturbation, is experimentally quite diffi- 
cult. This discrepancy raises the question about the per- 
tinence of the energy gain (ubiquitous parameter in tran- 
sient growth studies) as the most appropriate quantity to 
use for comparison between theory and experiments. 
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